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ABSTRACT 



A system is described which incorporated a Hewlett-Packard 
85 "desk-top" computer to control a frequency synthesizer 
and read the output of a lock-in analyzer :o measure, 
display and record the resonant frequencies, amplitudes, and 
quality factors of several modes of an acoustical resonator. 
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system is 


capable of locating. 


measuring , 


and tracking 


the 


resonan t 


modes as parameters 


which affect 


sound speed 


and 
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on are varied. An 


algorithm 


for rapidly 
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ing "good 


quality" measured 


data to 


a resonance 



line shape is described which determines quality factors to 
precisions cf better than 0.1 percent, amplitudes to better 
than 0.01 percent and center frequencies to better than 
O.lppm. Sample output is presided for the lowest three 
plane wave modes of a n air filled cylindrical resonator in 
the temperature range of -15 to 25 degrees Celsius. 
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I 



INTRODUCTION 



A. BACKGROUND 

Many applications in physios and engineering require the 
ability to determine the center frequency (fo) , maximum 
amplitude (A), and quality factor (Q) of a resonance. A few 
examples of this need arise in che measurement of the speed 
of sound, in the reciprocity calibration of elect roacou Stic 
transducers [Ref. 1], and in the measurement of sound 
absorption. 

Prior to the advent of computers ani instrumentation 
interfacing, severe limitations were imposed upon the 
experimental precision obtainable. The variation of 
amplitude at resonance is second order in frequency so that 
a small uncertainty in the maximum amplitude results in a 
greater uncertainty in the frequency corresponding to the 
maximum amplitude. Also, the Q is calculated from the ratio 
of a large number (center freqiency) to the difference of 
two large numbers which may differ by a small amount (the 3 
dB down frequencies), resulting in an uncertainty in Q which 
can be quite large. 
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This project was designed to produce an automated, high 
precision system capable of obtaining and tracking -he 
center frequency, maximum amplitude and 2 of several 
acoustical resonant modes in any given resonator as some 
external parameter affecting these values is varied. The 
approach was broken up into four task groups. 

1. Aco usti ca l Signal Generation and Detection 

If was necessary to determine which method of 
resonance excitation and signal detection was best suited 
for automated data acquisition and analysis. The following 
are the three techniques examined experimentally: 

a. swept frequency continuous wave excitation with 
tracking narrowband neterodyr.e reception using 
an HP 353DA spectrum analyzer. 

b. phase-sensitive detection utilizing a two-ohase 
lock-in analyzer and a programmable frequency 
synthesizer as a signal source. 

c. broadband noise excitation in conjunction with 
an FFT analyzer. 

2. Data Analysis Program Develo pmen t 

It was necessary to determine a precise and 
efficient computer algorithm for the fitting of the 
resonance data to a resonance Lineshape to extract fo. A, 
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and Q 



3. System I mpieme nt ation 

This involved the assembly of the required 
components (computer, analog-to- digital converters, counter, 
etc.) in a custom rack, interfacing and interconnection of 
components to the computer, resonator, and each other, and 
the writing of the software to allow the co mputer/contro Her 
tc orchestrate all of the acquisition, analysis and display 
functions. 

4. Svstem Test and Evaluation 

The final task was to test the system on the lowest 
frequency plane wave modes of a gas filled, cylindrical 
resonator tc determine: the precision with which amplitude, 

Q and center frequency could be measured; the tracking 
capabilities; the processing speed; and the utility of the 
output. 

B. CONCLUSIONS FROM TASK COMPLETION 

Upon completing tasks 1 and 2, it was determined that 
method b (frequency synthesizer/phase-sensitive detector) 
provided the best noise rejection and resolution, was 
computer controllable, and resulted in a data set which 
simplified the implen e neat ion of the resonance fit algorithm 
chosen from task 2. Method a of task 1 was rejected due to 
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the lack of a suitable 



computer interface which would allow 



automatic changes in frequency range. Also, the sweep w 
require catching the values "on the fly" thus resultin 
inaccuracies which were sweep rate dependent and data po 
which were not necessarily equally spaced. Method c was 
utilized due to the inability :o concentrate the energ 
the bands of interest, resulting in the need for hi 
power to the resonator without an increase in the sigr.a 
noise ratio in the frequency domain of interest. 

The completion of all of these tasks resulted in 
system to be described in the remainder of this thesis, 
system is capable of a precision equal to or better than 
part in 13 7 in frequency, one part in 13,000 in amplit 
and one part in 1000 for 2. 
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II* 12UI£^SE2 SETUP AND DESCRIPTION 

A. EQUIPMENT SETUP AND INTERFATE 

Figure 1 shows a o lock diagram of the equipment utilized 
in the final system. The following is a brief look, at the 
overall fun ct ion of the system which will be followed by a 
description of each individual piece of equipment.. 

The HP-85 computer directly controls all of the 
equipment, with the exception of the Lock-in Analyzer, via 
the Hewlett Packard Interface 3us (HPIB) . The computer 

sends a value for frequency and amplitude to the HP3325A 
Synthesizer /Function Generator, which causes an excitation 

of the acoustic resonances in tie resonator. The amplitude 
analog signal received by the Lock-in analyzer is converted 
to a D.T. voltage which is sent to the HP3497A Data 
Acquisition /Con trol Unit where it is digitized prior to 
being sent to the computer via the HPIB. After processing, 
the computer sends output to the HP2673A Printer and the 
HP7470A Graphics Plotter. The HP3456A reads the temperature 
in degrees Celsius from a thermistor mounted in the body of 
the resonator and, wtien interrogated by the HP-85, sends the 
current value via the HPIB. In the test system (Chapter IV) 
the HP3325A outputs the required frequency and amplitude 



DIGITAL VOLTMETER 
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FIGURE 1 BLOCK DIAGRAM OF SYSTEM 



sigr.ai by means of a power amplifier to the eiectret 
transducer at one end of the resonator. An elec tret 

transducer at the other end of the resonator picks up the 
signal and sends it via an amplifier to the E38G 5204. The 
ZG&G 5204 is phase-locked to tie HP3325A and outputs a D.C. 
voltage which is proportional to the Pythagorian sum of the 
amplified in-phase aid quadrature components of the signal. 
This analog amplitude is sent to the HP3497A which, when 
interrogated, sends via the HPI3, a digital measurement of 
the voltage level. The HP-35 also outputs values to the 
HP2673A and HP7470A to be printed and plotted respectively. 

B. EQUIPMENT DESCRIPTION 
1. HP- 85 Comjouter 

The Hewlett-Packard Personal Computer (HP-85) is an 
8 bit microcomputer containing as standard a 16K byte 
memory, a built-in thermal printer and a small CRT. The 
computer is programmed in BASIC computer language. As 
utilized in this application, the memory was expanded to 32K 
bytes through the use of an add-on 16 k byte memory module 
(#52903A). In addition, the HPIB interface (#82937A) , and 
the ROM Drawer (#62 936A) with a plotter/printer ROM 
(#00085-150 02) and ir.put/outpit ROM (#00085- 15 003) were 
utilized. 

1 5 



The HPIB is He wlet t- Paccar d ' s implementation of the 
IEEE Standard 488-1975 and is a parallel bus of 16 active 
signal lines, 3 data control lines and 5 interface 
management lines. It is capable of interconnecting up to 15 
instruments. The input/output and printer/plotter ROMs 
further expanded the memory in order to handle the 
transformations necessary to transmit information via the 
HPIB to the peripherals. 

2. Di g it a 1 Voltmeter HP3455A 

Although the HP3456A has numerous capabilities, in 
this application it was only utilized as a means to obtain a 
temperature from a thermistor and to transmit that value via 
the HPIB to the HP-85 computer. Its only utilization was to 
demonstrate the system's abilioy to track resonances under 
conditions of varying temperature. 

3. Syn th e si z er ^Fu no tion Generator HP3325A 

The HP3325A has a frequency range from 1 microhertz 
to 20 megahertz witi a resolution of up to eleven digits. 
The peak to peak output amplitude can be set from one 
millivolt to ten volts. It is fully programmable via the 
HPIB and was so utilized in this application. The execution 
time for a frequency command via the bus is 7.0 msec plus 



2.8 msec for each digit plus 2.3 msec for a decimal plus 
12.5 msec for each delimiter. The execution time for 
amplitude is 6.8 msec plus 2.3 msec for each digit plus 2.8 
msec for a decimal plus 13.0 msec for a delimiter. 

4. Data Acquisition^! ontrol Unit HP3497A 

The HP3497A was utilized as a voltmeter in this 
application. It received an analog voltage level from the 
lock-in analyzer which was proportional to the full scale 
meter deflection and converted it into HPI3 compatible 
signals to be transmitted to tie HP-85. It was completely 
program controlled, sampling and sending only when so 
directed by the HP- 35. 

5. EG8 G Lock-in Analyzer Model 5294 

Because the 5204 was not HPI3 compatible, desired 
functions had to oa manually controlled vice program 
controlled. This included changes to full-scale sensitivity 
and time constants. It received a phase reference 
synchronization signal from the HP3325A and the data signal 
from the resonator. Through vector manipulation it outputs 
the square root of the sum of the squares of the inphase and 
quadrature components of the amplified signal to the 
HP3497A. 
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5. Power Amplifier HP467A ini Sianal Amplifier. H?4 6 5 A 
These were ordinary amplifiers inserted in the 
system to provide signal gain. The HP457A was set for a 
voltage gain of two to drive tie resonator while “he HP465A 
provided a gain of +2) d3 for tie output from the resonator. 

7. Intelligent itaohics Printer HP26 7 3 A 

The HP2673A prints bi-airectionally an the rare of 
120 characters per second wirh a 7 X 11 dor matrix character 
font in a 9 X 15 character celL wirh twelve character sets 
to choose from and three modes of formatting. This printer 
was utilized to obtain standard sized paper output vice the 
small sized thermal printer output of the HP-35. It was 
utilized in the fully program controlled mods. 

3- Gra ph lcs Plotter 3? 7470 A 

The HP7470A was utilized to produce the high 
resolution and excellent quality graphs included in this 
report and suitable for publication. It is a two pen 
plotter capable of utilizing up to ten colors through 
programmed stops. It has a ciaracter plotting speed of up 
to six characters per second in any of five character sets 
with text written In any direction/ with or without 
character slant, and in varying sizes. It was utilized in 
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the fully programmed node via the HPIB. One of its aicest 
features is the ability to define the platting area, thus 
giving the operator tie capability of enlarging or shrinking 
the output size as desired. AIL figures in this thesis were 
produced on this piece of equipment. It is to be noted that 
although these figures are actual system outputs or replicas 
of system outputs, they were collected over a period of time 
which involved several experiments. Thus, each individual 
variable may not exactly coincide in value from figure to 
figure since parameters such as temperature, drive 
amplitude, full scale sensitivity, time constants, frequency 
band, etc., ail affect each variable and were possibly 
different for the figures shown. 

9. Thermistor 



The thermistor 


utilized in 


this set up 


1 5 


an 


HP-0837-0164 resistive 


type with a 


range of -80 


deg 


rees 


Celsius to 150 


degree s 


Celsius with 


an accur ac-y of 


plu 


s or 


minus 6 degrees 


Celsiu s 


for the range 


of -75 degrees 


Celsius 


to 130 degrees 


Calsiu s. 


Because th 


e temperature 


was 


only 


being utilized 


as an 


indicator of 


change, it 


WqS 


not 


necessary to 


measure 


this quantity 


accurately for 


this 


experiment . 


Thus, 


the most 


convenient me 


ans 


of 
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accomplish! r.g the desired taste w as utilized. The thermistor 
has a nominal resistance of 5303 ohms at 25 degrees Celsius 
rising to 3.684 Mohms at -80 degrees Celsius and falling to 
S2.7 ohms at 150 degrees Celsius. 

10. 8 e so n a tor 

A detailed schematic representation of the resonator 
utilized to test tae system is shown in Figure 2 [Ref. 2]. 
Although this system can be utilized with any resonator, a 
description of the one actually used will be given. 

The resonator is a cylindrical brass cavity with 
diameter and length egual to 2.54 cm. The ends are capped 
by simple electret transducers. A small slot (0.16cm X 
1.6cm) at the cylinder midplans exists to allow gas to enter 
and leave the resonator (a helium recovery line) as utilized 
in another application [Ref. 2]. 

A permanently polarized 12 micron thick disk of 
teflon (the electret), aluminized on one side, is 
electrostatically held against a sand blasted metal 
electrode, which forms the electrically active element in 
the resonator, while the aluminized side is connected to 
ground . 
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FIGURE 2 CROSS SECTIONAL VIEW OF RESONATOR 



The use of pre-charged electret material eliminates 
the need for bias supplies ordinarily present in 
conventional capacitive transducers. The colored-in circles 
in the figure are o-rings utilized in assembling the 
resonator to provide a gas tight seal. 
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A simplified flow chart of the program is shown in 
Figure 3. A brief overview of the flow will be given 
followed by a description of the individual blocks. 



Th 


e operator 


inputs an 


upper and 


lower frequency, a 


ve 


amplitude 


and a time 


con stant . 


This time constant 



will be discussed later, bur i: also needs :o be see on the 
iock-ir analyser as well as entered into the computer. This 
then constitutes the input to ’Search'. In 'Search', an 
approximate center frequency ana amplitude is displayed for 
all resonances within the frequency band entered. The 
desired modes are cm sen by tne operator and carried into 
•Sort'. In ’Sort' each of the modes is remeasured and then 
refined via a five point fit routine which results in an 
output of center frequency, maximum amplitude, quality 
factor, and bandwidth. A new time constant is determined 
and all are carried into 'Measure'. In 'Measure', a mors 
accurate center frequency, maxinum amplitude and Q are found 
for each mode. These then go into 'Ravine* where they are 
refined to a high precision by a modified Chi-squared 
minimization algorithm. The Last step is 'Track* where each 
of the modes is continually checked as well as the resonator 
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FIGURE 3 FLOW CHART OF PROGRAM 



temperature until a change in any resonant amplitude of 10 
percent or a difference in temperature of 0.2 degrees 
Celsius is detected at which time the program loops back 
into 'Measure' and repeats the routine. 

It is appropriate, at this time, to state that 
throughout the program, any system equipment response times 
(as mentioned in Chapter II) necessary for data 
transfer /acquisition are allowed to elapse through the use 
of WAIT statements prior to sampling any given response. 



A. INPUT 

The operator selects a start frequency and an end 
frequency to def in a the frequency spectrum to be 
investigated. Also, a drive amplitude and maximum drive 
amplitude are selected, which can be any value up to and 
including 3500 millivolts (limit of the HP3325A 
Synthesizer /Function Generator). Two amplitudes are 
selected because later in the program the drive amplitude 
will automatically increment or decrement the value in order 
to keep the lock-in output in its midrange. The second 
value (max. drive) is a default value adopted should the 
program attempt to exceed it while incrementing the drive 
amplitude. Finally, a time constant (tau) is selected. The 
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lcck-in analyzer divides the spectrum into 256 regions and 



integrates over each. The program calculates a maximum tau 
by 

r = 256/4 * ( f U pp e r -i OW er^ 

in order that there is nor a space between regions and chat 
the sensitivity over the region be nearly flat. This value 
is a maximum such that the operator selects a value 
available on the lock-in analyzer which is equal to or less 
than this and enters the value set on the lock-in, into the 
computer. 



E. SEARCH 

The frequency range is divided into 256 parts as 
mentioned above. Each of these frequencies is sent via the 
HFIE from the HP-85 to the HP3325A and a resonator amplitude 
value for each of these frequencies is obtained by the HP-85 
from the HP3497A via the HPIB. There is a programmed delay 
(WAIT) of 4 tau between the sending of a frequency and the 
sampling of its amplitude in order to allow the resonator to 
attain at least 98 percent of its steady state value as 
determined by the lock-in integration time. The choice of 
256 regions was dictated by the horizontal resolution of the 
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HP-85's CP.T display. These 255 values are then plotted on 
the CRT with the option of a hardcopy printout from the 
HP 7 4 7 0 A as in Figure 4. This graphically displays where 
each of the resonance modes within the specified frequency 
interval is located. 

Frcm Figure 4, it can be seen that not only are the 
responses from three major plane wave modes displayed but 
also chose of several other (azimuthal, radial, mixed) 
modes. If is possible to concentrate cn any mods displayed. 
After the operator enters a threshold relative amplitude 
value, the computer provides a printout on its thermal 
printer listing all amplitudes and their associated 
frequencies which equal or exceed this threshold value (see 
Table 2 in Chapter IV for printout). If the operator is 
satisfied with the output, then the plot and the thermal 
printout are utilized to input to 'Sort' the approximate 
center frequency and bandwidth for each desired mode. If 
the operator is not satisfied, the program loops back to 
'Input', new parameters are chosen, and the process repeats 
itself . 



29 



FIGURE 4 
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PLOT OF RELATIVE AMPLITUDE VS FREQUENCY (FULL SPECTRUM) 



C. SORT 



The purpose of 'Son' is 
the modes to be investigated 
values of center frequency, ma 
factor for each desired mode, 
is defined as sufficiently accu 
values to be automatically coat 
achieved by the operator input 
bandwidth for each mode desir 
configured, up to nine modes ca 
limit can be increased Through 
whatever number is desired. N 
convenience to simplify memory 

In this portion, the bandw 
into one hundred (an arbitrary 
resolution) equal frequency s 
sent and their respective ampl 
manner as in search. 

Once again, a plot is ge 
option for a hard copy (Figure 
satisfied with the output, ch 
the start of 'Sort* and new i 



to identify for the computer 
and to obtain "reasonable" 
ximum amplitude, and qualify 
'•Reasonable 11 , in this case, 
rate to allow these obtained 
rolled in ’ ileasurs*. This is 
ting a center frequency and 
ed to be worked upon. As 
n be selected. However, this 
a program/dimension change to 
ine was chosen as a matter of 
n anagemeni . 

iith of each mode is divided 
number which gave sufficient 
taps which are individually 
itudes collected in the same 

aerated on the CRT with the 
5) . If the operator is not 
ea the program loops back to 
aputs are entered. If the 



30 



FIGURE 5 
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PLOT OF RELATIVE AMPLITUDE VS FREQUENCY. MODE 



operator is satisfied, the program finds the maximum 
amplitude and its corresponding frequency (center 
frequency) . This is accomplished by examining the value of 
the amplitude of each of the one hundred points and 
selecting the one which is greatest. The amplitude value of 
the 3 dB down points is then found from this maximum 
amplitude. The program again sorts through all of the points 



and finds 


the closest 


"bra cket 


la g" 


po i nt s 


to this 


calculated 


3 dB down 


p cint . 


k li 


near int 


er po 


la t ion 


is then 


utilized to 


calculate 


the 


cor res 


ponding 


f r e 


quency 


for the 


calculated 



3 dB down amplitude point. Df oourse this process occurs 
twice since there is both an upper frequency and a lower 



frequency for t 


he 3 d3 


down ampli 


tude value. This process 


constitutes the 


Hf V 3 


point fit”. 




From these 


values. 


a quality 


factor is calculated which 



is equal to the center frequency divided by the difference 
cf the upper and lower 3 dB down frequencies. This value, 
as well as the center frequency and maximum amplitude, is 
output to the 3P-35 thermal printer (see Table 2 of 
Chapter IV for printout). 

This procedure is repeated for each of the remaining 
modes. Figures 6 and 7 ac e the plotter outputs of the other 
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PLOT OF RELATIVE AMPLITUDE VS FREQUENCY. MODE 



FIGURE 7 
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PLOT OF RELATIVE AMPLITUDE VS FREQUENCY. MODE 



two normal modes investigated in this example. Figure 8 is 
an example of one possible reason • for desiring no request a 
new input (wrong renter frequency). In this case, it is 
possible that one of t he 3 d3 down points could not be found 
from the data (cut of range). From this graph, the operator 
would obviously either choose a new center frequency input 
cr increase the bandwidth such than the 3 15 down point 

would occur in the data. 

D. NEW TAU 

In this block, three things occur. First, a new time 
constant is found. Phis time constant is the acoustic time 

constant required for each point to attain 63 percent of its 
steady state value. This is in analogy with the time 
constant required for a capacitor in an electrical circuit 
to charge up to its full value. The acoustic time constant 



for each 


mode is found by dividing 


the Q 


for the mode b 


y pi 


trines tne 


center frequency of th 


e mode. 




dcCi us 0 or 


the 


manner in 


which Q is found, this 


equates 


t 0 


dividing on 


e by 


PI times 


the difference of the 3 


aB down 


upper and 1 


over 



frequencies. After calculating the time constant required 
for each mode, the new system time constant selected is the 
greatest one required, thus insuring that the worst case is 
adequately covered. 
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PLOT OF RELATIVE AMPLITUDE VS FREQUENCY. MODE 



Second, -he background noise level is determined by the 
following method. The synthesizer outputs a frequency off 
resonance into the resonator and the lock-in analyzer is 
sampled one hundred times. The mean and the standard 



deviation of this set 


of measurements 


is calculated. 


The 


standard deviation is 


equated to the 


noise of the sy 


stem. 


It is utilized later 


in determining 


the signal to 


n oise 


(SNR) ratio for each 


mods. If is 


also utilized i 


n the 



calculation of the third item of this block which will be 
described next. 



The final 


(third) item is 


the examination 


of each 


mode 


determin e t 


he maximum number 


of points w hich 


car. be t 


aken 


h c u t h a v in g 


the possibility 


of the amplitude 


value of 


two 



adjacent points being equal or in the worst, case inverting 
position due to associated uncertainties or random system 
noise. If the frequency step size ( A f I is large enough, 
then 



A A - A ma v — A ( f o+Af ) 



can be chosen so 
Assuming a Rayleigh 



that the inversion 
distribution lineshape 



will 

then 



not 



occur. 
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A/Q 



A ( f ) = 




+ 



1 

Q 2 




Replacing f with (fc ► Af> and - = Af/fo, the resulting 
equation is 



AA « A 1 



l 

l(2€Q) 2 + i ] 1 ' 2 



For the system descrioed, 24 poinrs are normally taken in a 
frequency range of 2 5* rfhere 3 = fo/Bw. Therefore, ~hs 
following equation results: 



2€Q = 2 Af/3W = 1/6 < 1 



and, that quantity squared is <<1. Expanding and dividing 
by A 



A A 

T 



2 (€Q) 2 = 
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If A is chosen to be some number of standard deviations, a , 
based upon the noise measurement, then 



N <J 
A 



N 

SNR 




where °/A is 1/SNR (s ignal-to-noise ratio) . Now, if is 
found to be 



A f > 



1 



N 



2 ( SNR) 



Q 



If I is the number of points no be oaken in 2 3W then 



2 BW 2 f Q 

I 4 = — 

A f Q Af 

If the probability of an inversion, as described above, is 
set at 0.1 percent, then N = 3.3 based upon a standard 

distribution curve. It is worth noting at this point, that 

the probability of inversion away from maximum is everywhere 



8 ( SNR) 
N 
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lever until the 



~ h S 



order of 



signal level is on 



the 



background noise. utilizing the 0.1 percent criteria, 
I < = 1 . 55 * SQR (SNR) . This number (I) is output :o the HP-85 
thermal printer. The program then defaults no twenty-four 
points as a minimum. This value was based upon a 
calculation from the Nyguist sampling rate. The minimum 
from the Nyguist sampling rate (in integer value) would be 
three; but it is generally agreed upon in engineering 
applications that a factor of three or four is required in 
order to obtain a satisfactory reproduction of the signal. 
At this point the bandwidth was doubled to insure inclusion 
of the 3 <53 down points and thus the minimum number of 
points obtained equalled 24. For this choice the 

requirement that an inversion die to noise has a probability 
of 0.1 percent or less requires the SNR to equal or exceed 
237 (+47.5 <3B) . 



The system has worked wnen the indicated number of 
points was as low as 1 0 (still defaulting to a choice of 24) 
but occasionally gave faulty values due tc the reasons 
discussed above. Although it did not always happen, it is 
recommended that whenever the indicated number of points 
falls below 24, either the SNR is improved (increase drive 



amplitude, ezc.) or the mode is net suitable for further 
investigation. A sample of tie thermal output format for 
tau, noise and numoer of points, is shown in Table 2, 
Chapter IV. 

2. M2 AS UR 2 

2ach resonance mole has a frequency band calculated for 
it in ’Measure'. This band is twice the band defined by the 
3 d3 down bandwidta. Also, the relative value of the 
amplitude for the center frequency is evaluated to ascertain 
the need for adjustment in the drive amplitude. This is 
achieved by driving the resonator at the center frequency of 
each mode and sampling the amplitude value. If it is 
between 0.3 and 0.95 as measired on the Data Acquisitor. 
System, then no adjustment is necessary. If it is not 
within these limits, it is either doubled cr halved 
depending upon whether it was above cr below the limits. 
This new value is again tested in the above manner until it 
is either within the limits or has exceeded the maximum 
drive amplitude which was set back in 'Input'. If this 

happens (exceeding) , then the system defaults to the set 
value of maximum drive amplitude. This adjustment is 
necessary because the full scale sensitivity setting on the 



lock-in analyzer is a 
function . 

In ’Measure’ the 
the same manner (frequ 
’Sort', only this ci 
hundred are utilized 
allows the system to r 
value to fully utilize 
xhe H? 43 9 7A (as ex? 
TAU’ ) . Once the dat 
point fit is applied, 
closes x adjacent pcin 
define a parabolic cu 
parabola ana the corre 
setting the first der 
to frequency to zero, 
ana subsequent insert! 
equation yields a u 
frequency. This a 
calculation of the 
corresponding to the 
linear interpolation 



manual rather than computer conxrolled 



resonator amplitudes are obtained in 
ency sen "--amplitude measured) as in 
me, twenty-four frequencies vice one 
and xhe 'WAIT' becomes 12 iau which 
each within 6 ppm of the steady sxaxe 
xhe 5 1/2 digit A-to-D conversion of 



lain s d 


or - vioiisly at the end 


of 


' NEW 


a are 


collected, the follow! 


ng s 


even 


The 


maximum amplitude point 


a r.a 


the 


is o n 


either siie of it are 


use 


d to 



rve. The maximum amplitude of this 
spon ding center frequ ency are found by 
ivaxive of the amplitude with respect 
Solution of the resulting equations 
on of the solution inxo the original 
nique maximum amplitude and center 



mpli 


tude is then util 


izsd 


in the 


3 dB 


down amplitudes. 


7? 


requencies 


3d3 


down points are 


then 


found by 


in 


the same manner 


as 


pr evio usly 



4 2 



described in 'Sort'. A new quality factor is calculated 
utilizing these new fr eg uencie s and all three of these 
values; Q , maximum amplitude and center frequency are sent 
into * 3a vine'. 

?. RAVINS 

The values for matimum amplitude, center frequency and 
quality factor determined by the seven point fit described 
in the previous section ('Measure') are substantially more 
precise than those possible using conventional techniques 
employing a wave analyzer or spectrum analyzer under manual 
control. However, the high speed numerical processing 

ability of the digital computer allows a greatly expanded 
treatment of the data for only a minimal investment in 
actual processing time. 

Implicit in the assignment of A (amplitude), fo 
(frequency), and Q (quality factor) as the three parameters 
which characterize a resonance line shape (as is done 
throughout this thesis) is the assumption that the 
amplitude, as a function of frequency, near resonance can be 
described by; 
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A ( fn) = A/Q [( fn/f o-f o/f a) 2 + (l/Q) 2 ] 1 '* 



where : A(fn) is the resonance lineshape 

A is the calculated maximum amplitude 
Q is the calculated quality factor 
fo is the calculated center frequency 
fn is the measured frequency (24 per mode) 

Using all twenty-four points, it is possible to obtain the 
sum of the squares of deviations of the data from the 
resonance lineshape defined by a particular choice of A, fo 
and Q as: 

24 

S 2 (A, f o ,Q) = y~jAr- - A (fn)] 2 

n=l 

where An is the measured amplitude of the n (th) data point 
and A(fn) is as above. 

Conventional data analysis then accepts the "best 
choice" for A, Q, ana fo as that which causes S 2 to be a 
minimum [Ref. 3]. This minimization for a linear function 
is a straight forward algebraic process. However, an 
analytic least square fit to a non-linear function requires 
matrix inversions. This, in turn, not only requires large 
blocks of memory and computer time, but also can easily lead 
to errors due tc internal round-off or truncation of digits 
within the computer [Ref. 4], For these reasons, the 
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following method was devised and utilized for each of the 



resonant modes being 


i nvest 


igatad. 3r 


iefly , S 2 


was computed 


for small 


variations 


in A, 


f o and Q 


as well 


as for the 


originally 


computed 


values 


for these 


terms fro 


m the seven 


point fit. 


These 5 2 


t er ns 


were then utilized in 


a parabolic 


fit in an a 


ttempt tc 


find t 


hose values 


of A , f o , 


and Q which 



would cause S 2 to be a minimum. Because the results of the 
seven point fit were already quite good, the procedure, 
briefly described above and expanded upon in the remainder 
of this section, converges very rapidly to an astonishingly 
high degree of precision. This precision will be analyzed 
and discussed in detail in Chapter V. 



For the remainder of this paper, the sum of the squares 
of deviations (S 2 ) will be referred to as the C-G value, 
which was chosen for historical reasons. 

One at a time, each of the three parameters (A, fo, Q) 
in the resonance Lineshape (A(fn)) ace varied. The 
following is an example to clarify this process. 

First the C-G term for the calculated 2 is found. Mow Q 
is varied by +0.5 percent ana a second C-G term is found. 
Finally, 2 is varied by -0.5 percent and a third C-G term is 
found. These three C-G terms are used to define a parabolic 



t hs 



curve (Figure 9) and a new value for Q which minimizes 
C-G term is found. 

This new value of 2 repLaces that found in the seven 
point fir and becomes the new value no be inserted in the 
resonance lineshape. This time, however, A is the variable 
to be adjusted. Just as occurred with 2, C-G terms are 
found for A, A + D.23A, and A - 0.2%A. Again, these C-G 

terms are utilized in a paraoolic curve minimization and a 
new value for A is obtained which now replaces the A found 
in the seven point fit. Finally -hese new values of Q ana A 
are utilized to determine a new fo in the same manner. This 
completes one iteration of the process and thus provides 
these new values of 2, fo, aid A for the next iteration. 
Each iteration refines the precision of the values of Q, fo, 
and A and provides this new set of values for the next 
iteration. This process is repeated until the next 

iteration contributes a negligibly small variation in the 
C-G value [ Ref. 3], 

In the actual program utilized, the oerms were varied in 
exactly the order described because in the relative 
precision of the terms, Q was tie least precisely known term 
(historically, this is always true). Since this was the 
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FIGURE 9 RAVINE PARABOLIC CURVE 



case, varying it first resulted in the greatest effect. The 
next lowest precision term was A, so iz was obviously the 
next choice to be varied. Finally, fo, which was actually 
well known, was varied. Because of the relative precision 
in Q, A, and fo the variations were arbitrarily chosen at 
the indicated values (0.5i» of 2, 0.21 cf A, and 0.0051 of 
B») , as a matter of experience since these values reflected 
a sufficient variation for the ' ?. A V I N 2 ' technique to be 
efftctive. The better the approximations for Q, A, and fo 
put info ’RAVINS’, the smaller these values can be made. 

Open completion of the final iteration (program 
dependent) the three newest values for 2, fo, and A are 
utilized to find a final C-G term which is then printed on 



the 


system 


printer aid 


stored o.i 


a data tape. This 


C — j 


term 


is 


related 


to Chi-sq 


uared and 


may therefore be 


user 


:o 



determine the "goodness of fit" [Ref. 3], At this time, 
other parameters and variables are also printed on the 
system printer and stored on the data tape (rime, 
temperature, pressure, mode number, fo. A, Q, SNR, drive 
amplitude, sequential data tape item number, and C-G Ravine 
value) . These will be discussed at the end 
and a sample output will be shown (Table 1) . 



d8 



of Chapter IV 



If the r e a re 



then a 



additional modes to be evaluated, 
loop is made bask to ' Measure • and the next mode is 
evaluated in the sane manner as described above. Upon 
completion cf the final mode to be evaluated, the program 
shifts info 'Track'. 



G. TRACK 

In this portion the currently measured amplitude of each 
mode is compared with the measured value obtained during 
'Measure'. If the value differs by +/- 10 percent 
(indicating a change in some external parameter) , then the 
program loops back to 'Measure' to find the new Q, A, and fo 
for each mode. Also, the temperature of the resonator is 
measured and if a difference of 0.20 degrees Celsius is 
found then the system returns to 'Measure'. If neither the 
amplitude nor the temperature has met the oriteria necessary 
to return the system to 'Measure', then the 'Track' section 
continually loops through itself until either cf these 
criteria is met. The system remains in this endless loop of 
'Measure', 'Ravine', and 'Track' until the operator causes 
the system to halt, the data overflows the available 

storage space on the data tape (850 mode entries) or a 
catastrophic system failure occurs (loss of power, etc.). 
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The preferable method of stoppage is to utilize a 'PAUSE' 
command from the computer keyboard since this enables 
continuation or program restart as desired. It. is to be 

noted that any keyboard action will result in a halt. 

This, as well as any of the other mentioned halts (except 
•PAUSE' ) 
prior no 



may/will regains re-loading the program into memory 
re-running the system. 
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SAMPLE EXPERIMENT 



In order tc test the system’s ability to measure the 
stated quantities and to track the resonant modes as some 
external parameter was varied, an air filled resonator was 
cooled in a liquid nitrogen hath. This system was chosen 
for its simplicity and because the variation in sound speed 
wi oh temperature for an ideal gas is well known. The 
resonator was placed in a Dewar flask and a small quantity 
of liquid nitrogen was introduced no cool the rescnanor. 
Cnee the temperature staoilized, the computer system was 
activated and the lowest three plane-wave resonance modes 
were tracked as the system warmed. This tracking and 
recording of data proceeded unassisted by the operator for 
approximately twenty-two hours. At the end of this period, 
the resonator had reached room temperature and the process 
was halted. The printout was examined and the data tape 
was utilized to produce Figures 10 through 14 on the system 
plotter to demonstrate the utility of the system. 

Figure 10 shows a point by point plot of temperature in 
degrees Kelvin versus time in seconds. After time 10,000, 
the curve begins to appear as a solid line because the 
individual points become too closely grouped together to be 
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visually distinguished in the pLor. This likewise occurs at 
varying places in the other figures. Also, in Figure 10, 
there is a slight dip at rhe beginning of the curve because 
the system was activated as the temperature was still 
slightly decreasing vice stable as descrioed aocve. 

Figure 11 shows a plot of the log of the difference 
between a reference temperature (a: infinite rime) and each 
individually measured temperature, versus rime. This plors 
as a nearly straight line whose slope is a measure of r'ne 
exponential thermal eguilibrium time. 

Figures 12 , 13, and 14 show plots of the square of the 
normalized frequencies of eaci mode (fo/M#) 2 versus the 
absolute te mperarur e. For an ideal gas these plors should 
be straight lines. In a perfecr resonator, each successive 
resonant mode would oe an exact multiple of the first and 
the normalized frequencies would lie identically upon one 
another. In this experiment mods two is parallel to mode 
one but slightly higier in normalized value while mode three 
lies diagonally in-between mode one and mode two. 3ecause 
the exact physical properties of the resonator itself were 
nor being investigated, this was not persusd furrhsr. The 
slight curvature of the lines is probably due to the fact 
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FIGURE 10 
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PLOT OF TEMPERATURE VS TIME 



FIGURE 11 
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TIME (10*3 SECONDS) 



FIGURE 12 
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TEMPERATURE (K) 



that, the thermistor does not indicate absolute temperarure 
exactly. The slop? of the lines and the length of the 
resonator could be utilized zo determine the constants in 
the ideal gas equation of state. 

These figures (10 through 1 i l demonstrate a few of the 
possible means of displaying the data stored on the tape. 
It is intended than future users will determine what plots 
are needed and generate them accordingly. The program used 
to produce one of these plots from the data tape is included 
in Appendix 3. 

To enumerate what data is available on the tape, a 
partial reproduction (dot marrix outpur is non permitted in 
this Thesis) of the system printer output for another 
experiment is shown in Table 1. The rime shown is nhe 
average time at which the measurements were taken. The 

temperature is the average at which time nhe measurements 
were taken (in degrees Celsiusi . The pressure is sec to 
zero because a varying pressure was not to be evaluated, and 
thus, a means of deTecting pressure was not included in the 
equipment setup. The next entry in the line is the mode 
number for which the data applies, followed by the final 
ravined values for fo. A, and 2* The SNR is obtained from 
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TABLE 1 SAMPLE SYSTEM PRINTER OUTPUT 
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5587 23.842 0 3 20524.59757 1. 4381E-002 83.7679 17297 7. 3880E-006 



the division of A by -he value found for “he noise (as 
explained in Chapter III). The next value on a normal 
printout would be the drive amplitude of the system followed 
by the sequentially increasing number of the data tape item. 
These two items were dele red from this figure in order to 
conform to thesis requirements without requiring photo 
reduction. The whole line constitutes one data item so far 
as the tape limit of 350 data items are concerned. However, 
each item in the line is an actual separate entry within 
each of the data items. The final entry is the C-G ’Ravins 1 
value. 



Table 2 


is a reproduction of the cutpu 


t coral 


red 


from 


the HP-35 


thermal printer during the 


course 


of 


an 


experiment . 


For further amplification o 


f each 


of 


rhe 



entries, the reader is directed tc Chapter III, sections B, 
C, and D. The system printer was utilized for this figure 
again because dot matrix outputs are not permitted in this 
thesis. 
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TABLE 2 SAMPLE OUTPUT OF HP-85 THERMAL PRINTER 



AMPLITUDE IN VOLTS 


FREQUENCY 


.01082675 


6817.00 


. 00829425 


6896. 00 


. 01056175 


13769. 00 


. 00820725 


13848. CO 


. 00828725 


20405. 00 


.01466500 


20484. 00 


. 02054900 


20563.00 


. 01539075 


20642.00 


.00869400 


20721. 00 



MODE I CENTER FREQ IS 6830 AND AMP IS .0181355 Q IS 190.123117729 
MODE 2 CENTER FREQ IS 13756 AND AMP IS .0127195 Q IS 135.493980736 
MODE 3 CENTER FREQ IS 20540 AND AMP IS .02531475 Q IS 194.224445807 

TIME CONST > = a 86062488511 NEW TIME CONST IS 10 



The mean is . 000378561 The standard dev is 9. 797528548 4E-7 



NO. OF POINTS FOR MODE 1 IS 212 NO. CHOSEN IS 24 
NO. OF POINTS FOR MODE 2 IS 177 NO. CHOSEN IS 24 
NO. OF POINTS FOR MODE 3 IS 250 NO. CHOSEN IS 24 
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V. EVALUATION and recommendations 

This section describes the method utilized to determine 
the ’'goodness of fit" and the resulting precision. The 
program previously described was utilized in a steady state 



cordi tio n 


(no variation 


of 


e.c t ernal 


parameters) . 


A Iso, 


•Ravine* 


was programmed 


for 


f i f t e e n 


iterations per 


mode. 


with a pr 


intout of the 


^ / 


A, and 


fo values for 


each 



iteration. In this manner the change per iteration could be 
evaluated for each of the variaolss. As stated previously, 
•Ravine' is completed ( C-G minimized) when successive 
iterations produce a negligible variation in the value 
obtained [Ref. 3]. ? cr a given sample run, tne variation 

between the fourteenth and fifteenth iteration ((15-14) /IS) 
was approximately one in one miLlion (mode 1) to one in ten 
million (mode 2) for 3, one in two million (mode 1) to one 
in fifteen million (mode 3) for A, and one in six billion 
(mode 1) to one in one hundred thiry-seven billion (mods 2) . 
It is evident that the change per iteration is definitely 
small and thus the fifteenth iteration values were utilized 



as a base line for toe calculation of precision to follow 



TABLE 3 


COMPARISON OF 


PRECISION OF 


ITERATIONS 




MODE 1 


MODE 2 


MODE 3 




Q 


Q 


Q 


3rd/l5th 


1. 04E-3 


8. S4E-5 


1. 25E-4 


5th/ 15th 


3. 85E-4 


3. 31E-5 


4. 64E-5 




A 


A 


A 


3rd/l5th 


4. 66E-4 


2. 85E-5 


3. 99E-5 


5th/ 15th 


1. 72E-4 


1. C5E-5 


1. 47E-5 




f 


f 


F 


3rd/ 15th 


1. 2CE-7 


3. 62E-9 


l. 54E-8 


5th/l5th 


4. 45E-8 


1. 34E-9 


5. 71E-9 
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Each C-G ravins value for the third and fifth iteration 
was compared to the base line to determine the precision of 
the system. This as accomplished by talcing the ratio of 
the difference of the base line and the iteration to the 
base line value. Table 3 shows the approximate (rounded) 



values obt 


ainea for 


eanh 


by 


mode. 


In as — epnesenc 


s a 


best /worst 


precision 


of o 


ne 


in ten 


thousand/one in 


ons 


thousand f 


or Q, of 


one 


ia 


t hirtv 


thousand/one in 


tve 



thousand for J, and one in thrae hundred million/one in one 
million for fo for the third iteration. The program, as 
presented, stops after the third iteration because each full 
iteration requires 45 seconds on the HP-35 and that was the 
point chosen in tha precision/time trade-off for this 
particular application. This will be further discussed at 
the end of this chapter. The results for the fifth 



iteration (in 


th e 


same manner) 


clZ 


e one 


in thirty 


theusand/one in 


thre 


e thousand for 


Qt 


o a s in 


one hundred 


thousand/one in 


five 


thousand for A, 


and 


one in 


five hundred 



million/ one in twenty million for fo. 

Precision has been discussed since it is independent of 
the equipment utilized. Accuracy, on the ether hand. 
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involves net only the precision to which any given 
instrument measures as well as that of the overall system 
but also the accuracy of each measurement and the accuracy 
of the overall system. For example, this system will 

measure the given value of the temperature to within more 
than five significant figures oetween -90 degrees and 150 
degrees Celsius, yet, it is only accurate to within six 
degrees Celsius (a fuction of the thermistor). Thus, there 
is a dramatic difference between precision and accuracy. 

A. CONCLUSIONS AND a 3 COMB END ATI ON S 

The ability of this system to obtain 2 r A, and fo far 
surpasses that of any manual method. In addition, it allows 
the operator to perform additional functions as the 
experiment progresses, dees not make transcription/ 

interpretation errors as is ? o ssibl e/likely with manual 
recording of data, gives highly reproducible results 
independent of the skill of the operator, and allows for the 
resulting data to be to be easily presented in a variety of 
formats in a very short period of time. 

It was originally intended that this complete system be 
a prototype to then be transferred to a system utilizing the 
HP9836 rather than the HP-35. However, the equipment 
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arrived too late to be incorporated 



. One major advantage of 
the HP9836 is that it has an eigh t-to-or.e speed advantage 
over the HP-85. Speed is the reason that only three rather 
than five iterations of •Ravine' were utilized. Since each 
mode took approximately 135 seconds, 6.75 minutes elapsed 
from the start cf ' Measure', dole 1, to the next start of 
'Measure', Mode 1. Obviously, the greater the number of 
modes the greater the time required. This time diffference 
could be a major problem if the changing external parameter 
caused the resonance of the moles to shift significantly sc 
that the response no longer fitted within the chosen 
bandwidth. One method of overcoming this would be to 
incorporate a prediction of the next location based upon the 
change of the parameter. This was not incorporated in this 
program for two reasons. One, it was never necessary in the 
experiment run; and two, the shift to the HP9836 with its 
speed advantage would eliminate this lengthy time period {it 
could perform twenty-four modes in the same period it 
presently takes for three). 

Another future change couli be to work in 'Ravine' for 
the first set of data while taking the measurements for the 
second set cf data. This would be a better utilization of 
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computer time since there are currently periods when the 
computer is idle (WAIT periods and time constants) . 

Another improvement would be to replace the lock-in 
analyzer with a computer controllable model. This is, of 
course, a trade-off of price versus convenience/time since 
there is approximately a two-to-one ratio in equipment cost. 

There are also numerous computer programming techniques 
which could shorten portions of the program and/or 
associated run times. One obvious change would be to make 

CET presentations of graphs optional or, if the system under 
test was already well known, to make the option available to 
enter the program at the ’NEW IAU* point thus eliminating 



the 


need to 


establish 


the 


desi 


red 


mode data. This 


was 


not 


done 


in tne 


presented 


prog 


ram 


5 i n. o 


e it will be utiliz 


ed 


by a 


vari 


ety of 


users and a 


ach 


will 


be 


able/is encouraged 




make 



modifications for specific applications. 
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APPENDIX a 



REFERENCE GUIDE TO PROGRAM LISTING 



Lins 


Nos. 


Ds scri pt 


1 - 


251 


INPUT 


270 


- 640 


SEARCH 


650 


- 1690 


SO RT 


1 65 1 


- 1899 


NSW TAU 


1900 


- 26 7 6 


MEASURE 


2679 


- 4997 


RAVINS 


4958 


- 53 2 0 


TR ACK 
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1 DIM F 4(200) ,A (300) ,C( 300) 

2 CLEAR 

3 PRINTER IS 2 

4 OUTPUT 722 ;"F4M6” 

5 CISP " F.Ef-.C V E PROGRAM TAPE AND INSTALL DATA TAPE" £ E EE 
? 

6 CISP "WHEN COMPLETED , ENTEF I" 

7 INPUT Y 

8 IF Y=1 THEN 5 ELSE CLEAF 0 DEEP £ CCTC 2 

9 EKASETAFE 

10 CREATE "DATA" ,850 ,86 £ REWIND 

11 CLEAR y DEEP 1 intial sec up 

12 SETTIME 0 ,0 

21 CISP "Enter lower tree" 

31 INPUT FI i lcwer free 
41 CLEAR is EEEP 
51 CISP "Enter upper rreq" 
o 1 InPLT F2 i upper tree 

7 1 IF F 2< = r 1 THEN CLEAR £ BEEP £ GOTC 51 
SI T= 1/ ( 4* ( ( F2-F 1) /2 56 ) ) l time constant 
91 CLEAR £ BEEP 

101 DISP "Set time constant on 5204 at " ; T ; " or smaller" 
111 DISP "Wnen complete, enter value set (in i'illi-sec)" 
121 INPUT T 1 

131 IF T1>T* 10C0 THEN CLEAR £ BEEP £ GOTO lu 1 
141 CLEAR i? EEEP 

151 CISP "Enter amp cr criv.jng tree in mv ( =< 3 50 C ) RMS" 
161 DISP " DEC INAL VALUES ARE NOT PERMITTED" 

171 INPUT A I amp of treq 

181 IF A> 2 50 0 THEN CLEAR £ BEEP £ GOTO 151 

182 FOR N=1 TC S 
lo 3 A £ ( N ) =A 

184 NEXT N 

1S1 CLEAR £ EEEP 

201 DISP "Enter largest value of arr.p ever aesirea {.=<250 
Omv RMS) " 

211 INPUT A1 I max future amp 
221 IF A 1> 3 50 0 THEN CLEAR £ EEEP £ GOTO 201 
2 20 I tn.ese are tne values to set up tr.e 2325A 
231 A$= v AL3 (A) 

2 32 A 2 $ = " AM " 

223 A 2 5= "MR" 

240 OUTPUT 717 ;A2$,A$,A3$ 

241 WAIT 200 

251 F 3=C EIL i (F 2-F 1 ) / 25 6 ) ! smallest int >= banowictn 

252 i SEARCH 

253 i SEARCH 

254 ! SEARCH 
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261 I next group senes free tc 3325A, anc gets air.p rrcir. 

34S7A 

26 S CLEAR 

270 DISP "I a it working in SEARCH" 

271 FOR jn=0 TO 256 
261 F4 (N) =F 1 +N*F 3 

2o2 ! F4(N) is tree sent 
29 0 F 4 $ = V AL $ ( F 4 ( N ) ) 

2 91 H 1 $ = " H 2 " 

292 H 2 $ = " F R " 

3 00 GelPUT 717 ; H2 0 , F 4 S ,H 1 $ 

310 v.AIT 4*T1+10C 

320 OUTPUT 70S vT3" 

.3 30 FR IER 70 5 ; A(N) 

3 5 u N E X T N 

351 ! t 5 ARC to APE SCALu FACTORS 

352 SEEP 

3c 0 F 5=f I- . i* ( FZ-r j ) 

37 0 F o=F2-h. l x (F Z-F 1) 

360 GCLEAR t CLEAR 

3 £ 0 S CA LE F 5 , F 6 , — .1,1.2 
4 oO X A X I e G ,1000 , F 1 ,c2 

4 10 TAXIS F 1 , . 1 . 1 ,1.2 

411 FOR N=0 TC 256 

412 PE NLP 

4 1 3 PLOT F 4 ( N ) , A (Iv) 

414 NEXT N 

4 15 L I S t "IF YCL EcSIPt A COPY EinIER ] , ELSE i*..TEi : 2" 

SEEP 

416 INPUT Y 
17 IF Y = 2 THEN 477 
410 GCLEAR & CLEAR & EEEP 
4 IS LISP "ENTER FIGURE # TC EE PRINTED " 

420 INPUT LF 

421 PLOTTER IS 705 

422 PEN I 

4 z z C A Ll, t 1- .1* (r i-t 1) # F2-r5C0,-.. ^ , 1 . : 

424 XAXIS u,100U,Fl,F2 

425 YAXIS FI, .1,0, 1. 2 

426 PEN 2 

427 FOR N=u TC 250 

428 PEN UP 

42S PLOT F4 (N) ,A (N) 

420 NE XT N 

431 PEN UP 0 PEN 1 

4 32 LEIR G, SIN (SO.) § PENuP 

4 o o FOR X=F 1 TC F2 ST EP 2 000 & PENUP 

4 2 4 MO v E X , - . 1 3 

4 o 5 LABEL vALF(X) 

4 3b N U XT' X 0 PENUP 
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440 LClR U 

441 FOR 1=0 TO 1.2 STEF .1 t PENUP 

442 t.CvE Fi-.05*(F2-F1),Y 

44 1 L A c E L 'vALS(l) 

44 4 NEXT 1 @ FENCE 
4 50 LL I K C 0 PENCE 

4 51 AC V E F1+30G0,— .21 

45 2 Lf-bEL "FREQUENCY IE rib" 

4 53 LEIF 0 e PEN CP 

4 5 4 AC v E F j , - . 3 

4 b 5 EPS EL " PLOT Or RELATIVE APrLITOEL vS r EE^uhi'jCl (roLL 
SFECTRUN) " 

45b LEIF. 0, SIN (PC) C PEiilP 
4 5 7 AC v E F L - . G e * ( F 2 - 1 1 ) , . 3 
45o L^-l:EL "RELATIVE « F.P LI 'I U C i_, " 

45 y LEIF 0 a FEN LF 

46 0 PC VE f j +2000 , 1 . 2 
4b 1 LAEEL "FIGURE ";L5 
402 PSNLP 

4 7b PLCTTE F< IS 1 
4 / / v^-LiiAi*' ^ SEEP 

476 EISF "Cc icu cesire tc rerun exp vvitn cirrerent cara 
n.e ters" 

475 CISP "Ir sc enter 1, ctnerv/ise enter 2" 

4c0 INPUT Y 

450 GCLEAR e CLEAR F SEEP 
5o O If 1 = 1 THEN 11 

510 EISF "Enter cecisicn pcir.t rcr anpiituce (o iC 1.2)" 

520 INPUT A 4 ! CECISICN POINT 
5:0 CLEAR t utEP 

540 EISP "Enter ruil-scale sensitivity setting rrcn; tne 
5204 m vc Its" 

550 INPUT A5 ! SENSITIVITY SETTING 

56 0 PRINT "An.pl ituce in vclts";" Frequency" 

£70 PRINT 

5 o 0 FCR N = 0 TC 25o 
5£0 IF A (N ) <A4 THEN 6 30 
6 CO A (N ) =A5 *A ( N ) 

6 1C PRINT USING o20 ; A(N) ,F4(N) 

o20 IIuAGE IX ,£. ELL'DCCOD , 1CX ,EDDEDC .ED 

b30 NEXT N 

631 PRINT USING 632 

o32 I PAGE 3/ 

640 CLEAR @ BEEP 

641 ! SCRT 

642 I SORT 

643 ! SCRT 

6 50 EISP "Flcvy irany rr.cces ac ycu cesire tc track (FAX CF 

9 ) ■■ 
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ooC INPUT A ! NufcBEF OF NCCES 
b 7 b EOF N = 1 10 R 
bSG CLEAR fe cLEP 

690 D ISP "nnat is center tree tor ir-cce " ; N 
7 00 INPUT ME) ! CEE 'I ER FRET 
710 CLEAR is EEEP 

720 DISF "Wnat is rrecwiotr. rcr a.cce ";K 

7:0 DISP "EL'S T EE GREATER THAN " ;2*F3 ! F:2 IS EL 

7 40 INPUT C(E) ! FREwLIDTH 

7 5 0 EE NT N 

7o0 CLEAR 

761 ! Frcn. 77u tc 960 is tne nrst rcugr. try rcr ir.easuri 
no (SOFT) 

7 70 FCF L=i TC y 

77 L LISP "I air wcrKinc m SCR'I rcr it.cce ";L 

/ c 0 F 0 R N =0 ‘i G 1 u L 

7Sc F 4 (N ) =N ( L) -C ( L ) / 2+N *C (L) /ICO 

bCO r 4S=vALS(f 4 (t ) ) 

olo CUT PUT 717 ; 1 2$ , F4 9 ,n 1 $ 

S20 a AIT T2*Tl+lub 
o 30 OUTPUT 7b9 ;"vT:" 
o 3 1 a AIT 50 
b 4 0 ENTER 709 ; A(N) 

3 5 u NEXT K 
651 EEEP 

Cob GCLEAR is CLEAR 

b 7 0 SCALE E s ( b) - . 1*0 (L) , F 4 ( 1 0 b ) + . 1 *C ( L ) ,-.1,1.2 

cbO XAXIS b,C(L)/lb,t‘i(0) , r 4 ( 1 0 0 ) 

oSU YAXIS F** (0) , .1 ,0 ,1. 2 

90 0 c’GF N = 0 TC 1 0 0 

910 PE^LP 

S2C PL Cl f 4 (N) , A (N) 

930 NEXT N 

940 ROVE F 4 ( 0 ) ,-. 1 

941 LAEEL "RCCE ";L;" ”;F4(C);" TC ";F4(100) 

94 2 LISP "IF YOU RANT' A COPY ENTER 1, ELSE ENTER 2",& EE 
E? 

943 INPUT Y 

944 IF Y=2 THEN Sb7 

94 5 GCLEAF <s CLEAR is EEEP 

94b LISP "ENTER FIGURE if TC EE PRINTED" 

947 INPUT LS 

954 PLOTTER IS 705 

955 PEN 1 

956 SCALE F 4 ( 0 ).- . 2 *G ( L ) ,t4(lbG) + .l*0(L) ,-.3,1.3 
9 57 /.AXIS U ,C (L) / lu , F4 ( u.) ,F4 ( 100 ) 

9 55 YAXIS F4(b) , . 1 , 0 , 1 . 2 

959 PEN 2 

960 FCF N =0 TC 100 
9o 1 PENLP 
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vr vo vc vo vo vc vo vr. vr vo vr> vo vo vo vo vr ud 



Sfc2 PLOT F4 (N) ,A (N) 

So 3 NEXT N 

6 4 FEnLP fe PEN 1 
o 5 LEIF 0 , S IN ( SO ) 

6 o FOR X=f4(C) TO F 4 { 10 0 ) STEP G(L)/iC t PENLP 
o 7 l\C V E X , - . 1 : 

68 LALEL V A L $ ( X ) 

6S NEXT X 

7 0 LDIR 0 § PEN O P 

71 FOR Y = t TC 1.2 STEP . ] O' PENLP 
7 2 i\C V £ r 4 { 0 ) - . 0 c *0 ( L ) , Y 

73 LABEL vALS(Y) 

74 NEXT 2 

/ 5 L i- I F 0 0 r u N c t 

7 o NO vE r 4 { 0 ) -b . 1 *C ( L ) ,-.21 

77 LASER " FREQUENCY IN HE" 

78 Ll IF 0 p PENuP 

7 P MOVE F 4 ( 0 ) - . 1 *C ( L ) , - . 3 

80 LAEEL "PLOT CF RELATIVE AXPLITLEE vS FRLceLN CY , ..’SEE 
" • r 

I AJ 

881 LEIF 0, SIN (SO) & PENLP 
S b 2 r C V E F 4 ( o ).- . 1 7*C (L ) , . 3 
S 83 LAoEL "RELATI VE AJ PLITLLE " 

P 8 4 LEIF 0 5 PENLP 
S c 5 i-.C v E F ( 0 ) + . 1 *C { L ) ,1.2 

ico LAEEL "FIGURE ";L6 <“ PENLP t P LETTER IS 1 
Sc7 CLEAR A E EEP 

8co ElSr "^c yes -want tc c.. ur.ee anycniry . Ir sc enter I, 
else 2 " 

PbS INPUT ie CLEAR e uLur- 
PSC IF 7=2 THEN lZCu 

luoC LISP "Cr.ar.^e sensitivity? Enter 1 ror yes, 2 rcr no 

if 

1010 In PUT Y& CLEAR f SEEP 
10 20 It 7 = 2 ThLN It 6 0 
1 0 30 CLEAR t BEEP 

1o 40 LISP "Enter nee sensitivity settinu" 

10 50 INPUT A 5 
lOoC CLEAR t DEEP 

lo7C DISP "Cnange rreceictn? Enter 1 rcr yes, 2 rcr nc" 

lObC INPUT Ye CLEAR e SEEP 

10SC IF \ — i'H Ei- 1 1 1 . 0 

1100 LISP "Enter nee rreceictn" 

1110 INPUT C { L ) 

1120 CLEAR (5 EEEP 

1130 DISP "Cnanc,e center rrec? Enter 1 for yes, 2 rcr nc 

11 

1140 INPUT Y£ CLEAF g BEEF 

1150 IF Y=2 THEN 1180 

1160 D1SF "Enter new center tree" 

73 



1 1 7 0 


I 0. ALT £• (L) 






lloO 


GCLE’AR CLEAR 




IliiC 


GC'IC 771 






lly 1 


! cren 1200 


to 16 bO is 


ter tree arc. tne 


S* 




12 Oo 


C L r. A F f CcLil/E 




1201 


LISP “I £ii 


w c r tc l r.c, 


C fi I j 


121o 


K , b 1 , E 4 , h 1 , 


£2 ,ri 3 , h 


4=0 


1 - 2 o 


il 2 # L : = 5 O 






12 30 


for' i. = 1 'lo- 


1 c c 




1^40 


ir A(W)<A(h 


) 1 H f.i\ 


1 ^ b L 


12 4 1 


I f A ( A ) -A ( h 


) IK 00 


1142 3 


1 2 *7 2 


t iv I rJ 1 A ( H ) 






’ -1 C. i 
X - U 


ct o (L) =A ( i\ j 






1 2 0 0 


r o ( L ) = F 4 ( i< ) 






1 <. 7 o 


ti 






1 2 b c 


1*4 E A *1 lm 






12S0 


A 7 (L)=Ag(L)/SvM2) 




1 3 0 o 


t'C F is — 0 IC 


h 




1310 


If A 7 ( L ) =A ( 


l) THEN 


14 lu 


1 32o 


If A 7 ( L ) <A ( 


i‘< ) l ti a 1 4 


1 27 L 


13 30 


IF A ( :\ ) < t j 


TEEN l45o 


134 0 


E i=A ( N ) 






1350 


h I -i 






1360 


GCTC 1450 






1 3 '/ C 


If A ( r,' ) >02 


T ri c i\ 1 


L U 


1 3c o 


E2 = A (i.) 






13 SO 


ri 4 - = i\ 






14 0 o 


bCPC 






1410 


El f L = A ( X ) 






1420 


Hi f ri 2 = W 






14 30 


F 7 =r 4 ( n2 ) 






1440 


GC'IC 1 4 7 O 






1450 


lm L X 1 X 






1451 


X= (A / ( L ) - £ 1 ) / ( £ 2,-h 


i) 


1^6 C 


r 7 =■ a * (f ^ (112 


) ~r 4 (Hi) ) + r ( 


14 70 


FOR N = H IC 


loo 




14 6 0 


IF A 7 (L) =A ( 


14) THEN 


15 oO 


14 So 


IF A 7 ( L ) > A ( 


N ) THEN 


15-0 


1500 


IF A (N ) >t 3 


THEN 15 


30 


151o 


c 3= A (14 ) 






152o 


H 3 = N 






1530 


GCTG 1620 






1540 


IF A (N ) <=3 4 


THEN 1 


57C 


1550 


£4= A (N) 






1 56 C 


H 4 = M 






1570 


GCTG lo2C 






1560 


£3 , £ 4 = A ( N ) 






1 5SC 


fi 3 , H 4 =N 






1 b 0 0 


F o = F 4 ( H 3 ) 







irst trv at 



a no w ret SORT frcce 



cen 

;L 
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161C GCT'C 16 40 
16 2 u NEXT N 

Is 2 1 X= ( A 7 ( L ) — £ 4 ) / ( B ; — £ 4 ) 

16 30 F8=- t X* (F4 (H4 ).-r 4 (H .;) ) ) +F 4 (B 4) 

1 fa 4 C <* (L) =f fa (E) / ( Fd-f 7 ) 

16 5u t-RIUT "MOLE " ; L ; 

1660 PRINT "CENTER FFEv IS " ; F 6 ( L ) ; " AND AB.P IS ";A6 IL) * 
A 5 

16/0 PRINT "o IS " ;G(D 
1671 CLEAR 

16 7 2 PRINT USING 16 7 3 
16 7 3 I RACE 5/ 

1 6 c 0 N EXT L 

IbSC CLEAR t G CLEAR 

169 1 DISP "I cn we r King on tuv.e ccr.stt.nt" 

169 S ! Frcir. 1700 tc 1740 tee largest Tire Constant ~s ic 
one 

1 7uC FOR L= 1 TO M 

1701 T (L) =o. (L ) / (FI*t 6 (L) ) 

1702 NEXT L 
1705 T-T(l) 

1710 FCR L=2 TC M 

1/11 IF T>T (L) THEN 1720 

1712 T=T (L) 

1 7 2 o NEaT L 

1750 BEET & CLEAR 

17 55 PRINT " T I R E CONST >= " ;T*luC0 

1760 LISP "Set ar.a enter new tire ccr.star.t (in it Uln-sec 

) " 

1770 DISP "RUST BE GREATER THAN " ; T ; " sec" 

1760 INPUT 2 

1790 T 1=1 

1791 CLEAR 

1 792 PRINT "NEW TIRE CCNST IS " ; T 1 

1793 PRINT USING 1794 

1794 I RAGE 4/ 

1795 ! MEASURE 
179o ! MEASURE 
1 i9 7 ! R.E A S o R £ 

1 7 's 6 ! R.E A. S U F £ 

1 799 DISP "I air. wcrKir.c cn stancarc cev , 4 cr cts, & var 
iance " 

LSCG ! Frcir. loGG to Is 24 is trie calc cr the ncise, sc. c 
ev . , anc tne number cr cts. useo . 

1601 G , v»=0 

16 0 2 F=(Fb (2) -Fo ( 1) )/4+Fb( 1) 

1803 F4$ = VALS(F. ) 

1604 OUTPUT 717 ; H2 $ , F 4 $ ,H1 9 
IcO 5 WAIT Tl*12 + iuC 
180o CLTPuT 709 ;"vTl" 
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1607 WAIT 1000 
1 6 0 c cOP. w = 1 '1C 1 o U 
1609 Els TER 70S ; A (2 ) 

1 6 1 o Vv = vv +A ( N ) 

1611 NEXT N 

1612 w=n/ 10 G 

1813 FOR N = 1 TC^lCo 
1614 S=(A(N)-w)"2 

1815 G=G+S 

1 6 1 o i\ E X T X 

lo 1 7 G =A 6 X £^R ( G / £ 'j ) 

1816 v. = Vv * A 2 

1815 r F I X T "Tne near, is ";w 

le20 FRIXT "Tne scancarc cev is ";G 

1 c 2 1 Jr F 1 1. 'I USING 1 8 i i. 

1 o 2 I i'Ai e E 4 / 

1825 FOR L= ] TC R 

18 30 I (L) =1.6 2*SvR (A 6 (L) *A.5/G) 

18 35 PRINT "L-C. Of PC I NTS FOR nCDE ";L; " IS " ; I (L) 

16 4u If I(L)<1CU THEN 1c 50 

1g45 I ( L ) — 1 0 0 

1550 If I(L)>50 THEN I860 

Lb55 I(L)=50 

1860 I (L) =1? ( I |L) ) 

1661 I ( D =24 

lct>4 PRINT "NC. CiiCSEN IS ";I(L) 

1g 65 PRINT 
16 co NEXT L 
1867 CLEAR t EEEr 
1660 PRINT USING 1861 



loOl 1 1M Gt 4/ 

18c 5 PRINTER IS 701,132 
lco7 DIM f$. [20] 

186c FS="C-G FAVINE" 



155C PRINT USING 1651 



"TILE" , "TEf P" 



REO ", "AMR ,"C" , "SNR" ."C-ANF" , " if FT" ,F$ 

1851 INAGE 4 A ,5a , 4A ,6 a ,5 A , 4X , 3A ,6 X , 7A , SX , 4 A , 1 
, 5 X , 5m ,4 X , 4A , oX ,1c A. 



1892 PRINT USING lcS3 ; "(SEC)"," I C ) " , "(KZ) 
" (n.v) " 



Hi. _ II 

7T f 



4X , 1A , 

" , " ( V r 



165c I i v A GE 5 A , 4 X , 5 A , a 5 X , 4 A , 1 C X , 6 A , c OX , 4 A , / 
1656 PRINTER IS 2 



1655 CLEAR 

1500 ! Frcir. 1500 tc 4580 is tne calculsticr. tor tne 
ure inducing center tree anc Q 



1SC1 15=1 

150 2 FCR L= 1 TC X 

1503 DISP "I AM IN LEASERS FCR NOCE ";L 
15Go Fl=f6(L)-Fb(L)/w(L) 

1907 F2=F6 (L)+Fb (L)/w(L) 



C. f 
>X , :A 

D " , 

H.GdS 



76 



lSUb U= (F2-F L)/I (L) 

ISOS v-{ F2-F 1) /2 
1 S L U h = 0 

IS LI If A o ( L ) > . 3 THEN lsL7 
LSL2 AS (L) =2* Ah (L) 

LSL 3 Ao (L ) =2*A 6 (L) 

LSL4 IF A S ( L ) < A L THEN 1SLL 

LSL5 AS (L) =AL 

LSL6 GCTC LS 2 5 

Lb L 7 IF Ao ( L) < . S5 THEN LS25 

LSLc Ao (L) =Ab (L) /2 

ISIS A 6 (L) =A6 ( L) / 2 

L S 2 0 GCTC LSL 7 

LS 25 A S =vAL v ( A S (L) ) 

LS26 CoiPOT 7L7 ;A2v,A$,A.3$ 

LS 27 v,AIT 2 0 0 

Ls2o ENTEF 722 ; T6 

LS2S FCF N-C TC I ( L ) 

I S 3 u r 4 (N ) = F o ( L) +0 * N — v 
LS 3 L F % $ -VA LS (F4 ( £ ) ) 

IS 3^ h L $ - " h 2. " 

L S 3 3 H2S = “ Fk" 

LS34 CLTPuT 7 L 7 ;H2v,E4S,Hiv 
LS 3 5 WAIT 12*TL + LCG 
IS 3o OUTPUT 70S ;‘'VT3" 

IS 37 ENTEF 7 0S ; A ( N ) 

LS 3b NEAT N 

Ie3S ENTEF 722 ; T 5 

LS4G FCF N = L TC I (L) 

LS4L IF A ( N ) < A ( H ) THEN lS6b 

LS42 IF A ( N ) =A ( H ) THEN 1S43 ELSE L S 50 

L S 4 3 PRINT A ( N ) 

LS 50 K=N 
LS60 NEXT N 
LS7G X=A IH ) 

L2tC T=(A(K + 1) -A(H-L) )/2 
1 9 S C 2 = ( A ( h - r L ) + A (Pi — 1 ) — 2 * X ) / 2 
2 Uotf F6 (L) =- (I'/ (2*2) ) 

2G L o A6 (L) = X + 2 * F b ( L ) + 2 * F o (L) "2 

2G2G Fo (L)=F4 (H) +F5 (L) * (F4 (ti) -F 4 (H-L) ) 

2 u :G 5 1,54 , H L , E 2 , ii 3 , H 4 = 0 

2 o 4 U 22,53 = 50 

223U A 7 (L) =Ab (L) /SvF (2) 

2300 FCF N=u TC K 
2 3 LC IF A 7 ( L ) =A ( N ) THEN 2*10 
2320 IF A 7 ( L ) <A ( N ) THEN 2 370 
2330 IF A(N)<Ei THEN 2450 
2 3 40 5 L =A ( N ) 

2 3 50 H 1 = N 
e ~ o U EjCT C ^.4bC 
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2 370 IF A(N)>D2 THEN 2400 
2 3bC £2=A|N) 

^ S C ti z = ii 

ii 4 L C G b. 1 0 <.450 



24K EL , E2= A ( N ) 

242u Hl,fi2=N 
2 4 30 F 7 = F 4 ( H 2 ) 

^440 oOTC <.4/0 
24 50 N E X T N 

24 5 1 X=(A7(L)'-b3)/(E2-El) 

246 0 6 7 = X * ( F 4 ( ri 2 ) -6 -* (hi) ) •+■ t ^ ( n I ) 
2470 6 C F. x\ -H TG I (L) 

24c b I? A7(L)=A(x;) THEN 25 60 
24 0 0 IF A 7 ( L ) >A ( a ) Tr.EN 2 5 4b 
2 5 0 0 IF A ( N ) > £ j 'I htc. 2 5 30 
251b b3 = A(.\) 

2 5 20 H 3 = N 

2 b 30 GCTC 2620 

2540 IF A (N ) <=£ 4 ThCN 2b 70 

2550 £4=A(N) 

<.«60 H 4 = A 



2570 

2580 

2550 

2o bb 
2 o 1 0 
cfczU 
2621 
2b 30 
2b 4 u 
2b 7 3 
2b 76 
26 7 7 
2678 

26 7 5 
2660 
2b d 5 
2ob5 
2690 

27 b0 

2710 

2711 

2712 
2720 
2 730 
2740 
2 7 50 
2760 
277G 
2775 



GCTC 2d 20 
E5 , E 4 = A ( ft ) 
hi ,H4 = N 
F 6 = F 4 ( n 3 ) 

GC.TC 2o4u 
NEXT N 

X=(A7 (L) -ts ) / (t 3-E 4 ) 

F6=- ( X* (64 ( H 4 ).- 6 4 (tid) ) ) +F 4 ( xi 4 ) 

0 ( E ) = c b ( L ) / ( 6 a -6 7 ) 

PRINTER IS 701,132 
CLEAR 

1 
1 

! VARY C 
FOR E = 1 TC 3 
FCR K=1 TC 3 

DISP "I AW IN RAVINE FCR 0 OF WCDE 
0(1) =C(L) 

0 ( 2) =1 .005*0 (L) 

0 ( 3) =. 95 5 *0 (L) 

R=A6 (L) 

U=( 1/0 (X) ) ~2 
E (K) =0 

FOR N= 1 TC I ( L ) 

C (N) -R/ (J (K) *SCR( (F4(N)/F6(L) -F 6 (L 

D = (A (N) -C <K) ) ~2 

E (K) =E (K) +D 

NEXT N 

CLEAR 



)/F4 ( N ) ) * 2 + 0 ) ) 
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2 7 80 NEXT K 

2813 v = -( (E (2) -E ( 3) )/( 2* (£ ( 2) +E ( 3) ~2*E ( L) ) ) ) 

2614 Q (L ) =J ( 1)+V*.GU5*J (1) 

2676 ! 

2877 J 
26 7 6 ] 

28 7 9 J VARY ANP 

2680 CLEAR 

2681 FCR K= 1 TO 3 

266 S LISP "I AM IN RAvINE FCR AMPLITUDE OF i^GDE ”;L 
2 b SO u (1 ) =A6 (L) 

2900 J ( 2 ) = 1 . 0 t 2 *A6 ( L ) 

2910 0(3)= . 996*Ao (L) 

2911 U= (1A?(L) ) * 2 
2920 E ( K) =0 

2930 FOR N= 1 TO I(L) 

294 0 C (N) =0 (R) / (0 (L) *SCR ( (F4 (N)/Fo(L) -F6 (L)/F4 (N) ) ~2+U) ) 

2950 D = (A (N) -C (N) ) ~2 
2960 E (K) =E (K) +C 
2970 NEXT N 
2975 CLEAR 
2980 NEXT K 

3020 V = - ( (E ( 2) -E (3) ) / ( 2* (E ( 2) +£ ( 3) -2*E ( 1) ) ) ) 

3030 A.6 (L) =J f 1)+V* . 002*0 ( 1) 

317 6 i 
3177 ! 

1 1 / o 1 

3179 J VARY FREQ 

3180 CLEAR 

3181 FCR K= 1 TC 3 

3169 LISP "I AM IN RAvINE FCR FREQUENCY OF NCCE " ;L 
3190 a ( 1) =F6 (L) 

3200 0 ( 2 ) = . 005 *F 6 ( L ) /Q ( L ) + F 6 ( L ) 

3210 J( 3) =-( . U05*Fto (L)/Q (L) )+F6(L) 

3211 F.= A 6 ( L ) 

3212 U= (1/C (L) ) "2 
3220 E ( K ) =0 

32 30 FCR N = 1 TC I(L) 

324 0 C (N) =R/ (C iL) *SQR ( ( F 4 ( N ) /O (K).-O (K)/F4(N))~2+U) ) 

3250 D= (A (N) -C(K) ) ~2 
3260 E (K) =E (K)+D 
3270 NEXT N 
3275 CLEAR 
3280 NEXT K 

3 30 5 v = — ( (E (2).-E (3) ) / ( 2* (E ( 2) +E ( 3) -2*E ( 1) ) ) ) 

3 30b F6 (L) =0 ( 1) +\/* . 00 5* J ( 1) /<* (L) 

3307 NEXT E 
330b ! 

3 309 ! 
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3 310 I 

3311 ! LAST SLOT FA VINE * 

CL EAR 

3 314 LISP " I Ai'l IN LAST SHOT R " 

3 3 lo J =Ab (L ) 

3317 0= ( 1/w (L) ) "2 
3 3 1 S E = 0 

3310 FOR u=l TC I (L) 

3 32C C (N)=0/(w (L) *SWR ( ( F 4 ( N ) / F c ( L ) - F 6 ( L ) / F 4 ( N ) ) ~2 + 0) ) 
3 32 1 L= (A (N) -C (h ) ) "2 
: z 2 2 E —2 + L 
3:2 3 fct. XT N 

3326 0=b*A5 

3327 Y=J/G 

3 32 o D III G6 [ 1 32 J 
3 3 30 ASS ICt.4 1 TC “DATA " 

3 3oO 1 4= T I A £ 

3371 T7=(T5+Tb)/2 



: :cij t 3 = U 



3 .3 ^ L 


PRINT* 1 


,15 ; Th 


, T7 


,? 


5 ,L , 


6 0 (E) 


,0 , 


C ( L ) 


t Y i A L 


(L) ,15 , E 


34 Cu 


AbS 1 L; 1 ; 77 


1 TC * 


















3 4 u i 


b: v> = 


X , 3£ . 3C , 


2 X , 


E . 


3Lc , 


2X , 2C 


A • 

,-A 


i c 


32 . 52 


, 4 A ,2 . 4 L £ 


,bX , 


3 C • 4 Li , : X , 


: CC :b , P 


, 4L 


,5 


X , 3L 


i c a , L 


. 4C 


£ " 






3 0 2 


PFIN'I LS 


INC C$ ; 


T 4 


,T 


7 c r 


,L ,F 6 


(L) 


i 2 


(L) , Y 


,A 6 ( L ) ,15 



34 lb I 5=1 5 + 1 

342b E 4 $. = vAL$ (F 6 (L) ) 

~ 4 : t AS = v AL 6 (n c ( L ) ) 

34 4b CCIPCT 7 17 ;A2S,AS,A3$ 

3 4 5 C ,,AIT 2GC 

:4ob CbTPbT 717 ; H 2 S , F - 5 ,E 1 5 
347 b a A I J. 'I i*lb + lbb 
3460 CbTPbT 70 5 ;"VT3" 

34 50 ENT ER 70 5 ; Cc(L) 

4 575 PR IT. TER IS 2 
46cC CLEAR 

4561 NEXT L 
4962 PRINTER IS 701 
4S63 PRINT LSING 4564 
4 564 I P;A G E / 

458; PRINTER IS 2 
4^9 5 ! 

4956 I 
4557 1 

4596 ! TRACK 

50 00 T5=u 

5010 FCR L= 1 TC A 

£u 1 1 LISP "I API TRACKING POLE " ; L 
£020 F4p=vA LS (r 6 (L) ) 

50 3b AS=VALS (A 6 ( L ) ) 
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50 4 0 OUTPUT 717 ;A2S,Ae,A3$ 

50 4 1 a A II 2 C v 

5050 OUTPUT 717 ; H2 8 , E 4 8 ,h 1 8 
5Coo ir.AIT il* 16+ 111 
5Col CLTPu'i 7o9 ;"vi3" 

50 70 ENTER 7 l 9 ; A i 

5080 IF A9 <= . 9*C c ( L ) THEN 5110 

5090 IF A$>=1. l*Cte (L) THEN 5110 

5 lUU GCTC 52 4 5 

5110 T 9= 1 

5 1 2 u IF A 9 > . 3 1 n £ N 5 210 
3150 Ad (L) = 2*Ac (L) 

514 0 A £ = 2*A S 

5150 IF Ao(L)<Al THEN 5120 

5 1 C L A6 (L) =A 1 

5170 GCIC 5250 

52 lo IF A_9<.£5 THEN 5245 

5<i4.L Ac ( L ) —t \ c { L ) / c. 

52 5l A £=A£/ 2 

5240 GCTC 52 lo 

5245 CLEAR 

5250 NEXT L 

5 26 0 cNTER ; '18 

5270 IF AES (TS-T7) >= . 2 THEN 1 E 0 2 

5 2 S 0 IF T9 < 1 THEN SOIL 

550c PRINTER IS 2 

5 30 1 T4.-TINE 

o 3 lo PRINT "GCNE TO NEAScRE LcE It n ? PLI x ulE r i TINE 
4 

d32l GOTO l£C2 

89 S 0 CLEAR 0 GCLEAR 

6999 CISP "THE END" 

POOL END 
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APPENDIX 8 



A sample program is lists! for obtaining information 
from the data taps aid utilizing it in a figure or plot. 
Ths key spot in this program is az line 31, zhe reading of 
data from the tape. There are a maximum of 550 records cn 
ths tape. In this case, 15 is indicating which record i s to 
be accessed. Each of zhese records can be accessed in a 



random 


manner. Only 


ths actual value 


is stor 


si, net the 


specif 


ic variable z ha t 


represents it. 


Thus for 


example, if 


in the 


original szor 


age program zhe 


values F 


= 1, T = 3, 


G = 4, 


and B = 7 wars 


stored, only the 


actual 


values would 



be stored in the order given (1,3, 4, 7). when reading these 



values from the 




zhe 1st 


value mighz be 


called A, the 


second Z, etc. 


The important 


point is that i f 


the 4th entry 


in the record 


is das 


ired , f 


our values must 


be read into 


memory (read n 


values 


to get 


ath value) and 


whatever the 


reading program 


ca 11 s 


ths e n 


try, is what i 


t now must be 


referred to wit 


hin ths 


program 
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1 CLEAR 0 E EE F g C-CLEAR 

2 DISP "hCw PANY POINTS ARE BEING EVALUATED? ENTER NUMEE 
E" 

2 INPUT E 

4 CLEAE 

5 LISP "ENTER THE LOwEST T BMP DESIRED CK FLCT" 
o INPUT El 

7 CLEAR & EEEP 

o DISP "ENTER THE hIGHEST TEMP DESIRED CN PLOT" 

9 INPUT E 2 

10 CLEAE (§ BEEP 

17 DISP "ENTER THE MODE # DESIRED" 

16 INPUT G 
19 CLEAR 

21 DISP "ENTEE THE FIGURE 5 DESIRED TC EE PRINTED" 

22 INPUT E9 
2 2- CLEAR 

27 DIM T1 (651) ,.T2 (851) 

28 GCLEAE ^ CLEAR 

29 ASSIGN# 1 TC "DATA" 

30 FOP 15=1 TC E 

31 READ# 1,15 ; 14 ,T7 ,P 5 ,L ,F 6 
22 IF - L=G THEN 4 0 ELSE 50 

40 Tl( 15) = (F6/G ) * 2/ 10 "6 

41 T 2 ( I 5 ) =T7 + 2 7 3 . 16 

50 NEXT 15 

51 FLOTTEP IS 705 

52 PEN 1 

60 GCLEAE t CLEAR 

61 B7=FLCGR (fil+273. 16) 

62 E8=CEI L (E2 + 273. 16 ) 

70 SCALE E 7- 10 , £8+ 5 , 3 5 ,5 1 
60 XAXIS 40 ,10 ,E7 ,B£ 

50 YAXIS E7 ,2 ,40 ,50 

51 LDIR 0 , S IN (90) g PEN UP 

92 FCR X=E 7 TC E8 STEP 1C g PENUP 

53 MOVE X ,3 7. 5 

94 LAE EL VAL$(X) 

95 NEXT X 

96 LDIR 0 

97 FCR Y=4G TC 50 STEP 2 g PENUP 
9 8 MCVE B7-4 ,Y 

99 LABEL VAL$(Y) 

100 NEXT Y 

101 LDIR 0 g PENUP 

102 MCVE E7+10 ,36 

103 LAE EL "TEMPERATURE (K)" 

104 LDIR 0, SIN (SO) g PENUP 

105 MOVE E7-8 ,37 

106 LAE EL " (FREQ/M#) "2 IN KHZ~2" 
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107 LDIB 0 0 PEN UP 
lu8 MOVE £7-6,35 

1C9 LAEEL "PLOT CF FREW "2 VS TEiVP .PlGDE ";G 

110 LD IF; C 0 FEN UP 

111 MOVE 37+5,50 

112 LAEEL "FIGURE ";ES 

149 PEN 2 

150 FOR N = 0 TO 3/ 5-1 

151 P.=G+ 3* N 

152 PEN UP 

155 PLOT T2 (P) ,T1 (P) 

15+ NEXT N 
155 PENLF 

160 CLEAR is GCLEAF 
Ifcl REMIND 

162 LISP "ALL COl\E, I IRE TO CC ON" 

170 EEC 
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